Response to Zhang et al.

T wo recent reports by Zhang et al published in Nov (2011) and our group (Anand & Kai, 2012; published online in Dec 2011) have suggested a role for the conserved Tudor domain protein encoded by qin as a component of piRNA pathway. Both studies agree that qin plays a role in the piRNA pathway and is required to repress reteroelements but differ in some aspects of their conclusions about qin function. We reported that kumo, hereafter designated qin, homozygous females exhibited mislocalization of many piRNA pathway components from the nuage, a perinuclear structure where piRNA pathway components are concentrated in the ovary. We also reported a reduction in piRNAs targeting transposons (Supplementary Table 3 in Anand & Kai, 2012). In contrast, Zhang et al found no change in localization of piRNA pathway components to the nuage and observed an increase in Aubergine(Aub)Aub homotypic ping-pong in the qin/Df alleleic combination, which led Zhang et al to propose that qin is required for maintenance of heterotypic, Aub-Argonaute3 (Ago3) ping-pong (Zhang et al, 2011). In the accompanying Correspondence, Zhang et al challenge some of our findings. The key issue is whether the localization of piRNA pathway proteins to the nuage is compromised in absence of qin. The two groups used different alleles in the original reports. Zhang et al used qin and qin (transposon insertion alleles) in their original article (Zhang et al, 2011). In our hands, these alleles over deficiency did not behave as genetic null (i.e. females were not completely sterile right after eclosing and still expressed partial qin transcript). The fact that Zhang et al did not observe nuage mislocalization in those alleles may reflect hypomorphic status of those alleles. Our study was based on an allele, qin, which has a small deletion that removes a potential start codon of qin and fails to produce qin transcript. It behaves like a genetic null in trans to a deficiency (fully female sterile). Zhang et al now report that they did not observe mislocalization of Aub, Ago3 and GFP-Vasa when qin was placed over Df(3R)Exel6180, and suggest that observed defects in nuage localization could reflect a background mutation on the qin chromosome. Although not reported in our previous report, our result had been confirmed in ovaries from qin/Df(3R)Exel6180 females and Fig 1 shows the result of an independently repeated experiment. In stage 5-6 egg chambers of the heterozygous control, Tudor domain proteins, Tejas (Tej) and Krimper (Krimp), and the Piwi-family proteins, Aub and Ago3, are observed as perinuclear foci. By contrast, in qin/Df (3R)Exel6180 ovaries, the prominent localization of these proteins in perinuclear foci was perturbed. Krimp and Tej foci made larger aggregates in the cytoplasm (with more than 80% of egg chambers showing clear mislocalization), while Aub and Ago3 showed larger aggregates slightly away from nuage with a similar penetrance (>70%). We observed the same defects in later stage egg chambers of qin/Df(3R)Exel6180 ovaries (data not shown). These results confirm our original finding that qin is required for robust localization of these piRNA patwhay components to perinuclear nuage (Anand & Kai, 2012). We also note that our earlier report showed that expression of a full-length qin transgene under control of the germline driver, nosGal4, rescued the sterility of the mutant (65% hatching rate compared to sibling control) and the defects in nuage localization of several piRNA pathway components in qin homozygous female (Fig 2 and Supplementary Table 1 in Anand & Kai, 2012), further supporting that the observed defects were mainly caused by the loss of qin function. In their correspondence, Zhang et al also claimed that piRNA levels between qin homozygous and qin/TM3 ovaries were similar in the small RNA libraries published by both groups (GSE34728, Anand & Kai, 2012; Zhang et al, 2014). Our libraries contain a large number of non-coding RNAs derived from ribosomal RNAs and internally transcribed spacers and our mutant library contained almost three times more endosiRNAs and flamenco-derived piRNAs than control libraries (Supplementary Table 1, this response). In addition to our previous normalization with noncoding RNA (reads mapping to rfam database), we now compared piRNA levels between the mutant and heterozygote after normalizing them separately with endo-siRNA (Malone et al, 2009; Handler et al, 2011), snoRNA-derived noncoding RNA (Taft et al, 2009), and repeatderived 21-nt small RNAs. All normalization methods bring flamenco-derived piRNAs to equivalent levels between mutants and controls: flamenco-derived piRNAs are not much affected by germline piRNA pathway mutants (Malone et al, 2009), thus making them a good parameter to examine the raw data normalization. The new analysis confirmed approximately 2.5-fold reduction in cluster-mapping piRNAs and 1.5–1.9-fold reduction in transposon-matching piRNAs in qin homozygous ovaries (Supplementary Table 1, this response). This reduction is comparable to what we reported